www.gusucode.com > Lucas-Kanade optical flow method for 3-D源码程序 > Lucas-Kanade optical flow method for 3-D源码程序/code/demo.m

    % testdata contains two subsequent 3-D images.
clear; clc;
load testdataLK3D.mat;
[ux,uy,uz]=LK3D(img1,img2);

% Enhance the quiver plot visually by downsizing vectors  
%   -f : downsizing factor
f=5;
x=ux(1:f:size(ux,1),1:f:size(ux,2),1:f:size(ux,3)); 
y=uy(1:f:size(ux,1),1:f:size(ux,2),1:f:size(ux,3)); 
z=uz(1:f:size(ux,1),1:f:size(ux,2),1:f:size(ux,3)); 

[X,Y,Z]=meshgrid(1:size(x,2),1:size(x,1),1:size(x,3));
quiver3(X,Y,Z,x,y,z); axis([1 size(x,2) 1 size(x,1) 1 size(x,3)]);